Spectrum and correlation functions of a quasi-one-dimensional quantum Ising model 
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For a model of weakly coupled quantum Ising chains we outline the phase diagram and establish 
that well below the transition line the system has a remarkably one dimensional spectrum. We 
study the dynamical magnetic susceptibility and find a very rich spectrum with several modes 
corresponding to a fundamental particle and its bound states. The approach is based on Bethe 
ansatz and the Random Phase Approximation applied to the interchain exchange. 
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The Ising model is one of the generic models for many 
quantum and classical systems. Besides extensive ap- 
plications to spin systems (here we refer the reader to 
0, 13, 01), the model could also be used to describe in- 
teracting electric dipoles (like in systems with orbital 
degrees of freedom 4j) or arrays of interacting Joseph- 
son junctions (see, for example, 5]). The low dimen- 
sions are the most interesting since here the model ex- 
hibits strongest correlations. In one dimension the quan- 
tum Ising model (see Eq.© below) is exactly solvable 
by means of Jordan- Wigner transformation which con- 
verts the spin Hamiltonian into a Hamiltonian of non- 
interacting fermions. In general the spectrum has a gap 
which is closed when the transverse magnetic field is 
equal to the exchange integral {g = 0). The exact solu- 
tion also exists when the magnetic field has z-component, 
although only when the x-component is kept equal to J 
0. The solution predicts a rich spectrum with as many 
as eight particles and a hidden Eg symmetry. An interest- 
ing question is whether some of this fascinating physics 
may survive in realistic systems which are almost never 
truly one-dimensional. The goal of this paper is to discuss 
this and related problems. We begin by introducing the 
model, go on to outline the phase diagram and then dis- 
cuss spin-spin correlation functions in the ordered phase. 
In certain regions in parameter space we see that the 
one-dimensional behavior is clearly visible. 

Let us consider the quasi-one-dimensional (weakly cou- 
pled one-dimensional chains) quantum Ising model de- 
scribed by the following Hamiltonian: 

H ^Y.^iD+ E ^™(*-jK("MM (1) 

= -J|iE{'^I(^K(" + i) + (i + 3K(")}j^2) 



where i,j label the chains, n,m label the sites on the 
chain and a are the Pauli spin operators. For simplicity 
we discuss the case where J,mi(i — j) = Jj_{i—j) if n = m 
and otherwise. 

The dimensionless parameter I + g is the strength of 
the transverse field written in a way such that the one- 
dimensional model has a phase transition at T = at g = 



h = 




QCP 



FIG. 1: Expected Phase Diagram of model Q. Estimates of 
Qc and Tc{g — 0) are given by Eg. 1141 and Eg . IIUI respectivelv. 
In the vicinity of <; = point the physics is well described by 
model 0. 



0. At r = the phase with 5 < is ordered with (cr^) / 
and the phase with g > is quantum disordered (see e.g. 
Ij). We would expect a non-zero inter-chain coupling to 
extend the ordered region to finite temperatures and shift 
the critical coupling to g 7^ as shown schematically in 
Fig. CI 

The spectrum and dynamics of the model close to the 
transition line is very well understood ,7| . The transition 
itself falls into the 3D Ising model universality class, the 
Quantum Critical Point (QCP) falls into the 4D Ising 
model universality class. We concentrate our attention 
on the region of phase diagram well below the transition 
line where new non-universal physics can be found. We 
demonstrate that non-trivial physical effects are possible 
in this region. 

If we are not very close to the phase transition line, 
the inter-chain coupling can be taken into account in the 
Random Phase Approximation (RPA) . The RPA expres- 
sion for the dynamical magnetic susceptibility in the dis- 
ordered phase is 

X(c., fc;ki) = [xi^iu;, k) - J±{k^)]-\ (3) 

where xiD is the susceptibility of a single chain: 

XiDico.k) = -iV rrfte""([a^(i,n),a-(0,0)]), 

(4) 
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where v = J\\ao is the on-chain velocity and oq is the 
lattice spacing in the chain direction. Formally RPA is 
the first term in an expansion in the inverse transverse 
coordination number l/z±, where z_l < 6 if we are em- 
bedding our chains within a three dimensional structure. 
The application of RPA to other quasi-one-dimensional 
magnets 8,'9] have demonstrated that far from the tran- 
sition it provides a good accuracy even for small coordi- 
nation numbers. 

In the ordered state expression 10} has to be modified. 
Namely, one has to replace xiD by the dynamical suscep- 
tibility calculated in the presence of an effective magnetic 
field generated by the neighboring chains other 
words, in calculating xiD one has to use the following 
Hamiltonian: 

Hid - E {--^11 [<^'in)a'i^ +!) + (!+ + ha'in)} 

n 

(5) 



real 

imaginary 



with the self-consistency relation 



(6) 



There are three energy scales in this problem. J|| the 
on chain coupling, m ^ \g\J\\ is the spectral gap for a 
single chain, and J± the inter-chain coupling. In order 
to treat the model as weakly coupled chains we must 
have J± <C J\\ and in order to apply the non-perturbative 
results for the individual chains which have all been found 
in the continuous limit, we must have m <C J||. However, 
m and J± can both be of the same order of magnitude. 

In the scaling limit the one-dimensional problem is in- 
tegrable along the lines g = or h = 0, and conformally 
invariant when both these conditions are met. Through- 
out this work we use the standard conformal field theory 
normalization of the magnetization operator, fixed by 



Xid{t,x) = {(T{T,x)a{0)) = 



(ao/v 



,1/4 



11/4 



^1-0, (7) 



where = + (x/u)^, r is the Matsubara time. 

We first calculate the phase diagram of the model 
plotted schematically in Fig. db). To estimate the crit- 
ical temperature, we set h = and then look for lo = 
divergences in the correlation function which signify a 
developing instability. 

We begin by considering the line g = where the un- 
coupled chains are critical. At finite temperature the 
spin-spin correlation function ^ becomes 



Xid{t, x) 



&mh{-KT{x/v — it)) sinh(7rr(a;/i; + it)) 



1/8 



(8) 



The Fourier transform Q gives the dynamic susceptibil- 
ity 113 

xMoj = 0,fc = 0) = ^{2T:Tao/v)-"^sm^B\l/lQJ/S) 

(9) 



FIG. 2: The three particle contribution to the dynamical sus- 
ceptibility for a single chain. 



where B{x, y) ~ T{x)V{y) /T{x + y) is the Beta function. 
Using the RPA we extract the transition temperature: 



TJJ« = 2.12 



z \ J \ 



Ji 



4/7 



(10) 



Now let us estimate the position of the Quantum Crit- 
ical point on the g axis. Again we use the RPA equation 
substituting in the expression for the dynamical spin 
susceptibility of the off-critical Ising model. According 

to [12 

XiD[r) ^ Ko[mr) + 0[e ) (11) 

TT 

where Zq — 1.8437. Hence for small uj, 



XiD{(^,k) = 



V Zo{mao/v)^^^ 



{vk)^ 



(12) 



where the dots indicate further terms in the form-factor 
expansion [l^ . The next term, the three particle contri- 
bution can be represented as|l4| 



Xi(s = 0) 
1 f dx f dy 

3 J Att J 47r (s/m)2 — [1 + 4 cosh a;(coshx -I- cosh?/)] 



f{2x)fiy + x)fiy-x) 



fix) = tanh2(a;/2). 



(13) 



with = — (vk)^ and is plotted in Fig. [3 The 
imaginary part of this term (which is the contribution to 
the structure factor) and all higher terms in the series 
is below lu — 3to, however all terms in the expansion 
contribute towards the real part at a; = 0. Fortunately 
the contribution from higher order terms are negligible 
at small s (see e.g. ^3)- We see explicitly in this case 
that xsis = 0)/xi(s = 0) = 0.0002 so that lO is indeed 
a very good approximation for small s. 

From Eq.El we can see that the QCP where Tc —> is 
given by the condition 



5c«1.42(ziJi/J||)4/^ 



(14) 
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In the vicinity of QCP the low-energy behavior of the 
quantum Ising model is universal and falls in the uni- 
versality class of the (d-|-2)-dimensional classical Ising 
model, where d is the number of transverse dimensions. 
Since d =2 corresponds to the upper critical dimension of 
that model, in three dimensions the fluctuations give only 
logarithmic corrections to RPA. The RPA expression for 
the susceptibility at the QCP is 



(15) 



where vj. = (1/2)^05^/^ J||dV_L(k_L = 0)/dki. 

Near the critical point, we can examine the shape of the 
phase boundary. For low temperature, the correlation 
functions at finite temperature can be well approximated 
Qby 



where is the quantum dephasing time given by 



1 = ^e-™/^. 



(16) 



(17) 



The RPA then gives the condition for a singularity at 
bj — 0, k — Q as 



llTl = ZoJ\\J^g^'\ 



(18) 



so the transition temperature in the vicinity of the critical 
point is approximately 



ln(l/y)-lnln(l/y) 
y = |(Zo(J^/J||).9-^/4-l) 



1/2 



(19) 



For 5 < 0, uncoupled chains are completely ordered at 
T = 0, and at T > and m not too small there is order 
on a length scale 



(20) 



A crude estimate for the transition temperature can be 
is given by 



which gives 



v/m, 



2m 



ln(TO/ Jj_) 



(21) 



(22) 



Now we go on to calculate the dispersion and spec- 
trum within the ordered phase, where we have to consider 
/i 7^ 0. Model (O is exactly solvable at g = only j^. As 



rrii/mi 




1.000 


0.247159 


1.618 


0.0690172 


1.989 


0.0209579 


2.405 


0.0122653 


2.956 


0.0021898 


3.218 


0.0011328 


3.891 


0.0001623 


4.783 


0.0000055 



TABLE L The masses and weights of the particles of the one- 
dimensional quantum Ising model in a magnetic field (after 
[31) The first three particles lie below the incoherent contin- 
uum which begins at tJ = 2mi . 



we have mentioned, the spectrum consists of eight par- 
ticles which for completeness are listed in Table |l| The 
contributions to the dynamical susceptibility for small uj 
decline quickly with the growth of the particle mass and 
the decrease in spectral weight of the mode. Therefore 
the magnetic susceptibility at T = can be well approx- 
imated by keeping only the first three poles below the 
incoherent continuum 



Xid{s) 



Ami 



2 3 

i=l 



mt 



(23) 



157rJ||/i^ 

with s defined as before. In this model, the mass scale is 
given bv[l7| 

mi(/i)/J|| = ai(VJ||)«/i5, 

ai = 4.40490858 (24) 

and the single-particle expectation value is 

{a^m = a2(VJ||)'/'^ 

a2 = 1.07496. (25) 

The self-consistency relations |51 therefore give 

h/J\\ = [a2J±(0)/J|]''/" 
m/J|| = ai [a2J±(0)/J||]*^^ (26) 
In the RPA |(2J) the dispersion is given by the condition 
J_ ^ Y ^ ' 

J| ^ s2 _ ^2 



Jl = 



Ami 



where the final relation is obtained from the self- 
consistent relations and 0^3 = (4/157r)^af /q;2 = 0.130. 
To second order in J± this is solved to give 



2/2 2 / 2 I /7 T 

s I m^ = TOj /mi + ZiJ± 



EZj mi J± 




(28) 



In numerical values this gives for the first three modes 



1.0 + 0.032Ji(ki)/J_i(0), 
{s/mf^ { 2.618 + 0.009Jj.(kj.)/J_L(0). 

3.956 + 0.003J_i(ki)/Ji(0). 



(29) 



The relative weights of each mode are not significantly 
changed from those for the pure one-dimensional case. 

It is remarkable how weak the dispersion in the per- 
pendicular direction is, so the ordered phase probably 
remains very one-dimensional in character even when 
J^/ J\\ is not very small. 

We can obtain some results slightly away from the in- 
tegrable lines. In the neighborhood of the QCP where h 
is small, the dynamic susceptibility is given by l(T^ with 
the mass replaced by [is) 



m->m[l + {h/J\\f +0{h% 



(30) 



For small h, we must also have 



(o-) = X(w = 0, /c = 0)/i (31) 
which combined with the mean-field condition Q gives 



h 

m/J\\ 



1 



\/.g[J||/(ZoJi(0))] 
{Z,J^{Q)/J.f'\ 



4/7 



- 1 



(32) 
(33) 



Here, the dispersion in the perpendicular direction is 
much stronger than in H29|l . It is given by H15(l at 
the QCP; in the vicinity of Qc the gap is given by 
^ i9c ~ g)^ in the ordered phase to be contrasted 
with m ^ [g - Qc) a± g > gc- 

Our results indicate that some of the beautiful physics 
of quantum Ising model in magnetic field with a hidden 
Eg symmetry may be observed even in a realistic quasi- 
one-dimensional model in its ordered phase far from the 
transition line (in the vicinity of point on Fig. 1). At 
this point it may be possible to observe at least three 
coherent peaks in the dynamical magnetic susceptibility 
whose relative strength is approximately 1 : 0.28 : 0.09 
(see Eas. ()23|l and Table 1). At this point the spectrum 
is remarkably one-dimensional in character. We expect 
that when one moves along the g-axis of phase diagram of 
Fig. 1 towards the QCP, the excitation gaps will decrease 
and the transverse dispersion will grow. At the QCP the 
spectrum is three-dimensional gapless and the spin-spin 
correlation function is given by Eq. H15(l with logarithmic 
corrections. These corrections will convert the pole into 
a continuum; at criticality all excitation modes will col- 
lapse into it. This is illustrated schematically in Fig. |3| 

A good candidate experimental compound to see these 
effects is NaV205. In jl3| it is shown that the charge 
ordering at low temperatures is well described by the 
Hamiltonian |2Jl where the a variables correspond to an 




FIG. 3: A schematic diagram showing what happens to the 
coherent modes as a function of g. To the left of the figure, 
we have three coherent modes with very one-dimensional dis- 
persions. As we move to the right, the transverse dispersion 
grows and the overall mass scale decreases until all modes 
collapse at criticality. 



isospin that describes which rung of a ladder the electron 
is sitting. These isospin excitations should be experimen- 
tally probable with optical absorbtion [2(il |. 
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